Numerical evidence of regularized correlations in spin foam gravity 
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We report on the numerical analysis of the area correlations in spin foam gravity on a single 4- 
simplex considered by Rovelli in PRL 97 (2006) 151301. We compare the asymptotics and confirm 
the inverse squared distance leading behaviour at large scales. This supports the recent advances on 
testing the semiclassical limit of the theory. Furthermore, we show that the microscopic discreteness 
of the theory dynamically suppresses and regularizes the correlations at the Planck scale. 



A fundamental theory of quantum gravity is expected 
to improve the UV behaviour of the non-renormalizable 
perturbative quantization of General Relativity. The lat- 
ter should nonetheless emerge in the low-energy limit, 
where it can be considered as an effective field theory. 
Among the key features to be reproduced is the dis- 
tance dependence of the free graviton propagator encod- 
ing Newton's law. In Loop Quantum Gravity (LQG) and 
its covariant version, the spin foam formalism [l| , the UV 
behaviour is expected to be cured by the discreteness of 
spacetime at the Planck scale. While the regularizing ef- 
fect of such discreteness is clear, it is rather non-trivial 
how the discreteness smooths out to a low energy limit 
given by the effective theory of gravitons. This is the 
problem of the semiclassical limit in LQG and spin foams. 
In this letter we give numerical evidence that a behaviour 
consistent with the graviton theory at low energies is dy- 
namically suppressed and regularized at the Planck scale. 

The study of the semiclassical limit has received a great 
deal of attention over the last few years and irnportant re- 
sults have been obtained 0. In particular, in 0, i, H, 0, 0| 
a program was started to compute the graviton propaga- 
tor using correlations between geometric quantities. An- 
alytic results show that the leading order at large scales 
is consistent with that of the free propagator from the 
linearized quantum theory, thus providing an important 
piece of evidence for the correctness of the limit. Here 
we confirm this result numerically, supporting the ap- 
proximations made to deal with the complexity of the 
calculations. Furthermore, the numerical approach also 
allows us to study the small scale structure of the prop- 
agator, where the non-perturbative effects of spin foams 
are dominant, and we show that the discrete structure 
does regularize the typically divergent behaviour of the 
correlations, suppressing them at the Planck scale. Fi- 
nally, we point out the limitations of the model used, and 
discuss the relevant developments to be made. 

We consider the area correlations for 4d Riemannian 
quantum gravity defined and studied in Q]. These are 
correlations between fluctuations of the areas around 
a given background g of a 4-simplex, and corre- 
spond to some components of the graviton propagator 
Gfj,vpaix,y) = {0\hf,„{x)hj,„{y)\0). For extensive motiva- 
tion and discussion, see [a] . We index the ten triangles of 
the 4-simplex by /, and the area eigenvalues are given by 



Ai = ip{2ji + 1), with jl half-integers. For simplicity, we 
choose the background to be a regular four-simplex, with 
all ten areas having the same value, Aq = £p{2jo + 1). 
Given two triangles a and b, we consider the following 
area correlation. 
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where J\f = J^j, '^giJi] is the normalisation, is 

the boundary state and K the propagation kernel, or 
path integral amplitude. These are model-dependent 
quantities, and we describe our choices below. The quan- 
tity h.{ja) = {A^ — Aq^/Aq represents an area fluctua- 
tion, or equivalently, the fluctuation of the metric ten- 
sor h^i, projected along the normal to the triangle. 
WabUo) is the spin foam discrete analogue of the pro- 
jections Gab = n-i^n'^n'^n^ Gfii,pa{a,b) of the continuum 
graviton propagator around a flat background, with the 
two points taken to be the centers of the triangles. If the 
theory has the right semiclassical limit, the leading or- 
der of ll]) should match the leading order of Gab, namely 
the free propagator corresponding to the linearized the- 
ory, which we recall scales as the inverse squared distance 
between a and b. 

The boundary state represents a dynamical co- 

herent state peaked around both the (canonically con- 
jugate) intrinsic and extrinsic 3-geometry of the classi- 
cal background q chosen. Given the discrete 4-simplex 
considered, this means choosing a background conflgu- 
ration for the areas and their conjugate variables, the 
dihedral angles. Taking the simple choice of the equi- 
lateral conflguration {jo, 9), all triangle areas are given 
by Aq = £p{2jo + 1) and all inside dihedral angles by 
= arccos(|). The explicit form of such a state is not 
known in the full theory, but control can be gained by go- 
ing to lowest order in the perturbative expansion. Based 
on analogies with the continuum linearized theory, in [J 
the following Gaussian ansatz was made, 

= cxp I - ^ ai,n Sji Sj,n + i9^{2ji + 1)|, 

(2) 

where aim is a 10 by 10 constant matrix and 5ji = ji — jo- 
The matrix aim is non-diagonal, but the symmetries of 



2 



the equilateral background reduce the number of inde- 
pendent entries to three (see discussion below). Having 
the squared width proportional to jo guarantees that in 
the large jo limit the state ^ is peaked around both 
conjugate variables of the background q = {jo,0) (e.g. 
Q). The physical boundary state for the full theory is 
expected to have a lowest order contribution, correspond- 
ing to the free theory, given by ^ with a definite aim- 
So if we know the full state, we can fix aim looking at 
its perturbative expansion. In this context, the parame- 
ter for the perturbative expansion can be taken to be jo'. 
as geometric areas are given by ipjo, the limit jo oo 
drives the £p expansion. Therefore, the parameter jo has 
a double role: on the one hand, it describes the back- 
ground geometry of the boundary; on the other hand, is 
the parameter of the perturbative expansion. 

The kernel depends on the spin foam model chosen. 
As in we consider the Barrett-Crane (BC) model 



^M-n(2j/+i)'{ioj}, 



(3) 



where the integer k parametrizes the choice of face weight 
in the measure, and the lOj-symbol {lOj} is an SO (4) 
invariant tensor constructed with Clebsch-Gordan coef- 
ficients. To study its perturbative expansion, recall that 
when all spins are homogeneously large, i.e. ji = Nki, 
N i—f oo, the lOj-symbol has a stationary phase contri- 
bution of the form 



fi{ji) COS {SR[ji] + (f>) , 



(4) 



where ^^{jl) is a non-oscillating function scaling like 
A^"'/', SR[ji] = Y.i{2ji + mUi) is the Regge action for 
a single 4-simplex with triangle areas Ai = £p{2ji -I- 1) as 
independent variables, and 6 is an irrelevant phase. It 
was shown in [llj (see also [l^l) that this is masked by 
a non-oscillating contribution D{ji) that scales like N~'^. 
This dominant contribution corresponds to a degenerate 
geometry for the 4-simplex. However, as conjectured in 
[J| and proved in Q, D{ji) is negligible in the evalua- 
tion of quantities like ([T]), where its non-oscillating na- 
ture fails to properly match the phase of the boundary 
state ([2]) . The work of this paper supports these analytic 
calculations, and thus is the first (indirect) numerical ev- 
idence of Therefore the asymptotic behaviour of ^ 
is effectively given by the kernel for Regge calculus with 
measure lli{2ji + 1)'= fi{ji). 

The emergence of Regge calculus in the semiclassi- 
cal limit of LQG and spin foams has often been advo- 
cated [isj , and would provide a solid bridge to low-energy 
physics: Regge calculus is a discrete representation of GR 
known to reproduce the correct linearized quantum the- 
ory thus if in ([T|) we use and K from linearized 
Regge calculus, we expect to obtain the right free gravi- 
ton propagator. This suggests a possible strategy to fix 
(Xim by evaluating the boundary state in linearized Regge 
calculus [Tsj . In general, the matrix aim is non-diagonal, 
and this makes evaluating ([1]) extremely challenging. To 



simplify the numerical analysis we study the case where 
aim = OL Sim is diagonal. Indeed, using the formula for 
the {10j}-symbol as an integral over S0(4), this choice 
allows us to perform each sum over ji separately, which 
simplifies the numerical task. The integral is then dealt 
with using Monte Carlo methods. We will come back to 
this point at the end, and discuss improvements in the 
boundary state. 

We then have a model with a single free parameter a, 
determining the width of the Gaussian. For different val- 
ues of a, we study the matrix ([T]) as a function of jo, in 
particular to test its asymptotic behaviour and support 
the analytic calculations which appeared in (3] , which we 
briefly recall here. First of all, fl^ {2ji + 1)'^ M(j;) is a mea- 
sure term in ([T]), thus it will affect only the higher order 
corrections ■ The leading order is obtained by approxi- 
mating ^ with (the second moment of) a Gaussian inte- 
gral, like in the continuum linearized theory, with action 
Qlm 5ji Sjm given by the quadratic term in the boundary 
state and the second derivatives of the Regge action. 
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The Hessian matrix of the Regge action on the equilat- 
eral configuration was computed in Q, and it is a 10 
by 10 matrix with a regular structure inherited from the 
regularity of the equilateral 4-simplex. Fixing a triangle 
a, there are only three distinguishable choices for b: the 
case when b is the same triangle, the six cases when b 
is an adjacent triangle sharing an edge with a, and the 
three cases when b is an opposite triangle sharing only 
a vertex with a. Correspondingly, for each row of the 
Hessian there are only three different entries. 
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where for each row fim ~ —9 occurs once, 7/2 six times, 
and —4 three times. We can then write 



Qi 



1 

jo 



Aim, Aim = (a Sim + i Hlr, 



(7) 



Notice that Him. is not invertible, just like the ori gina l 
continuum term, due to diffeomorphism invariance [15| . 
On the other hand. Aim is invertible for any a > 0, thus 
the boundary state effectively provides a gauge-fixing for 
the propagator. 

Following the same procedure as [3| , we obtain for the 
absolute value of the leading order 



|M^iiO(a)| = l|A-^,| 
Jo 



(8) 



If the theory has the right semiclassical limit, ([8]) should 
give the free graviton propagator of linearized quantum 
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gravity.^ Choosing the harmonic gauge, the latter is 
given by 

^ I -]\ 1 ^ lingua ~t~ ^lyp^utT ^niy^pa 

G,.p.(a,6) = -- ^^-^p . (9) 

For the equilateral background wc have chosen, there are 
only three independent projections Gab, and furthermore 
the squared distances are all proportional to jo- Both 
features are matched by ([5]), which has only the three 
independent entries discussed above and scales as 

In general, |Vl^^'^(a)| scales as (ajo)^^ for both a :s> I 
and a <^ 1. To have explicit values, we choose a and b to 
be opposite triangles, and fix a = 0.5 and a = 5. From 
(HI) we obtain 



|W^ip?(0.5)| = 



1.02 

jo 



0.13 

Jo 



(10) 



In Fig. [T] we compare these analytic results with the nu- 
merical analysis of the full formula ([T]) , where in the ker- 
nel (121) we choose the simple measure term with k = 0. 
The dots are the numerical evaluations of the absolute 
value of ([T|) for these cases, and the error bars are one 
a. Remarkably, good agreement is reached already at 
Jo 50. This is our first result: we numerically tested 
the inverse squared distance asymptotic behaviour of the 
area correlation in spin foam gravity. This is crucial to 
support the approximations used in [1] to deal with the 
complexity of the calculation, in particular the conjecture 
that the D term drops out. 

Let us now focus on the short scale behaviour. We see 
that the divergent behaviour of the graviton correlations 
gets regularized at high energies by the discrete structure 
of spin foams. This is how the full theory enhances the 
effective field theory where the latter breaks down, and it 
confirms the intuition that spacetime can not be consid- 
ered as fluctuating around the flat metric at the Planck 
scale. The peak is very close to the Planck scale, and 
its exact location depends on the value of a, but also on 
the measure term in (|4]). To show this, in Fig. [2] we plot 
again the case a = 5, but this time with fc = 1. Notice 
that the position of the peak is slightly pushed to the 
right, but the qualitative behaviour is the same. 

We want to conclude with an outlook for further devel- 
opments to turn this picture into a concrete prediction of 
the theory. This will require enhancing both the bound- 
ary state and the kernel used here. 



We take the absolute value for a better comparison with the 
linearized Gab, because the spin foam kernel (|4]| provides the 
complex exponential of the Regge action even in Riemannian 
signature. The phase of JHJ is irrelevant for our work. 
One of the three cases is the correlation of a triangle with itself. 
While the continuum graviton propagator between the same two 
points is divergent, this is not the case for Waa- This should 
not be surprising, as it is among the regularizing effects of the 
discrete microscopic structure. 
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FIG. 1: Numerical study of ^ (dots), versus the analytic 
result of the leading order, on a log-log plot. Top panel: the 
case Q — 0.5. Bottom panel: the case a = 5. Raw data and 



more plots are available at http://jdc.math.uwo.ca/graviton 
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FIG. 2: The case k — 1, a — 5. The leading order is the same 
as for A; = 0, and the position of the peak shifted to the right. 



The first improvement would be to remove the ambi- 
guity in aim by fixing it with a dynamical requirement. 
To obtain a quantitative matching with the free gravi- 
ton Gab wc need the right boundary state, with the right 
non-diagonal structure and values of the entries, coming 
from the dynamics and a choice of gauge. As mentioned 
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above, a possible way to obtain this quantity is to evalu- 
ate it in linearized area Regge calculus. 

Furthermore, from analogies with conventional quan- 
tum field theory we expect there to exist a procedure 
for extracting the full boundary state from the non- 
perturbative kernel. Extending this procedure in back- 
ground independent quantum gravity (see for instance 
jlGj ) would be extremely useful. The full state will fix 
ctim by looking at the perturbative expansion, but also 
contribute to the small scale structure, thus aS'ecting the 
numerical analysis presented here. This will certainly 
modify numerical results such as the exact location of 
the peak, but will most likely not change the qualitative 
picture of suppressed correlations at the Planck scale. 
This regularizing effect is very likely to survive in the full 
theory, because it is a facet of the discrete microscopic 
structure, more than of the details of the model. 

On the other hand, the kernel itself needs improve- 
ment. Indeed the BC kernel does not capture the dy- 
namics of quantum gravity in a fully consistent way. Here 
we used it in the restricted context of area correlations 
on a single 4-simplex, where it provides a sensible quan- 
tum gravity amplitude. To test the full tensorial struc- 
ture of ([5]), we need to consider also projections that in 
the discrete picture would correspond to correlations be- 
tween the dihedral angles of the boundary geometry. As 



pointed out in [5|, |17[ , these correlations can not be stud- 
ied using the BC vertex amplitude. Furthermore, the 
calculations presented here need to be extended to many 
4-simplices, and the large spin limit of BC ^ precisely 
lacks the constraints which are needed to correctly treat 
the areas as independent variables. These issues have re- 
cently been addressed, and promising new models have 
been proposed [Tsj . It would be extremely interesting to 
apply the same analysis to these new models. 

In conclusion, we have presented a spin foam model 
where correlations consistent with the graviton theory at 
low energies are dynamically suppressed and regularized 
at the Planck scale. This work shows that the spin foam 
framework for quantum gravity naturally provides the 
expected regularization of the corrections at short scales, 
i.e. high energies: the theory has a short length scale 
appearing dynamically, which also suggests that space- 
time can not be considered as fluctuating around the flat 
metric at the Planck scale. We stress that the model pre- 
sented here is not predictive, and further work is needed 
before the correctness of the semiclassical limit of spin 
foam gravity can be claimed. Our results give a glimpse 
of what the qualitative picture of quantum gravity and its 
bridging to low-energy physics could be like, and suggest 
interesting new questions for further investigation. 
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